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Abstract 

We investigate the dynamics of electron spin qubits in quantum dots. Measurement of the 
qubit state is reahzed by a charge current through the dot. The dynamics is described in the 
framework of the quantum trajectory approach, widely used in quantum optics, and we show 
that it can be applied successfully to problems in condensed matter physics. The relevant master 
equation dynamics is unravelled to simulate stochastic tunneling events of the current through the 
dot. Quantum trajectories are then used to extract the counting statistics of the current. We 
show how, in combination with an electron spin resonance (ESR) field, counting statistics can be 
employed for quantum state tomography of the qubit state. Further, it is shown how decoherence 
and relaxation time scales can be estimated with the help of counting statistics, in the time domain. 
Finally, we discuss a setup for single shot measurement of the qubit state without the need for 
spin-polarized leads. 
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I. INTRODUCTION 



Controlling and preserving coherent quantum dynamics in the framework of quantum 
information processing is a challenging task.— Very recently, more and more experiments 
on implementing such ideas in mesoscopic systems based on solid state devices^ have been 
realized, e.g. Josephson junctions,-^^ and also single electron spins in single defect centers.— 
The electron spin in quantum dots has been recognized early as a potential carrier of quan- 
tum information,^ but experimental developments of suitable mesoscopic devices have only 
recently been pursued. 

In previous work it was shown how quantum dots may serve as spin filters, or memory 
devices for electron spin.— Important progress was made in both theoretical and experimental 
research focusing on measurement schemes through charge currents . '^''^' ■ 'i^'^i^^i^^'^i^^i^^'^i^^ Even a 
single-shot readout of the electron spin state has been realized^ and allows for the measure- 
ment of the relaxation time of a single spin. Still, the decoherence time of a single electron 
spin in a quantum dot has not yet been determined experimentally. 

In some of these experiments important quantities are the counting statistics of tunneling 
electrons.— As for charge qubits, a measurement on the single electron level may be achieved 
through a single electron transistor (SET) device^^i^S or with a quantum point contact (QPC) 
close to the quantum dot .^^i^^ 

It is important to realize that the measurement through a charge current itself has dy- 
namical implications for the measured qubit. We are thus lead to the problem of noise 
and statistics induced by the measurement process in these mesoscopic systems.— Such 
problems have been tackled some time ago very elegantly through the concept of quantum 
trajectories in quantum optical applications.— In particular, jump processes to describe 
the time evolution of open systems while counting emitted quanta are well established in 
the framework of systems that are described by a master equation of Lindblad type. Such 
ideas have already been applied to measurement processes based on quantum point contacts 
in mesoscopic devices.—"^ In the context of quantum information processing, such quantum 
trajectory methods turn out to be essential for the design of active quantum error correcting 
codes^l and, more generally, of quantum feedback mechanisms.— 

Historically, a major driving force behind the development of quantum trajectory meth- 
ods were the growing possibilities to experiment with single quantum systems in traps or 
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cavities. More recently, such experiments have been extended to mesoscopic sohd state 
devices. Therefore, we expect a growing need for such methods in these fields. 

The aim of this paper is to show how quantum trajectories serve as a useful framework 
to discuss the physics of mesoscopic carriers of quantum information under continuous mea- 
surement. In particular, we determine counting statistics of electrons tunneling through a 
quantum dot, depending on the electron spin-state. We show how a simple setup for state 
tomography can be achieved through a measurement of counting statistics in combination 
with a coherent ESR field. We display how decoherence and relaxation time scales can be 
extracted from the measured data, in the time domain. 



II. ELECTRON SPIN DYNAMICS ON THE QUANTUM DOT 

We consider a quantum dot with spin-| ground state in the Coulomb blockade regime as 
inS*2iifi^ see also FigO The quantum dot is subject to a constant magnetic field Bz which 
leads to a Zeeman splitting = gfisBz of the electronic states, where g is the electron g 
factor and fiB the Bohr magneton (throughout this paper we use \g\ = 0.44 for GaAs and 
units such that h = 1). Two leads at chemical potentials fii and /i2 are coupled to the dot 



for charge transport. Further, as in Ref. |lO|, we allow for an electron spin resonance (ESR) 
field to drive coherent transition between the two spin states. 

Leaving sources of uncontrollable environmental infiuences aside for a moment (see be- 
low), the total Hamiltonian consists of contributions from electrons on the dot, electrons in 
the leads and a tunneling interaction between dot and leads, 

-f^tot = -f^dot + -f^icads + Ht- (1) 

Here, i^dot = Hq + -ffESR(^) contains contributions from charging and interaction energies 
of the electrons on the dot, the interaction energy —^A^cTz with the static magnetic field, 
and the ESR-Hamiltonian -ffESR(^) = —\gi^BB{t)ax of the interaction of the electron spin 
with a magnetic field B{t) = B^cos{ujt — ip), oscillating linearly in the x direction. The 
(Tj {i = X, y, z) denote the usual Pauli spin matrices. The Hamiltonian for the two leads 
{k = 1, 2) reads i^ieads = J2kna ^kn(^kna^knai with c^^^ tlic crcatiou Operator of an electron 
with orbital state n, spin a and energy eun in lead k. Finally, the coupling between dot and 
leads is described by the standard tunneling Hamiltonian i/x = J2knma TknCkncrdma+ h.c. , 
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FIG. 1: Closed dot (left): Chemical potentials are too small to allow an electron to tunnel onto 
the dot. Open dot (right): After excitation of the dot electron, the chemical potential /xi is large 
enough for an electron of lead one to tunnel onto the dot and form the singlet state with the dot 
electron. 

where we denote with a tunneling amplitude and with dma the annihilation operator of 
an electron on the dot in orbital state m. Following Ref. for the description of the dot 
dynamics in Sect. Ill Al we will also include further (microscopically unspecified) dissipative 
interactions between the dot states and their environment that are not among the known 
contributions to the total energy as they appear in Eq. ((H). 

In the following we give a qualitative picture of the relevant dot states, see FigHJ more 
details may be found in Ref. [lOl For simplicity we assume there is only one electron on the 
dot. With cr^lt) = +11)5 the electron has the ground state |t) with energy = Q and the 
exited state ||) with = A^. If an electron tunnels onto the dot, the two electrons can 
form the sing let state \S) = (Hi) - |iT))/V2 with energy Es or either of the three triplet 
states. As the triplet state |T+) = |tt) has higher energy (for suitable magnetic field^^), 
the singlet 15) is the ground state for two electrons on the dot. The chemical potentials 
are chosen such that Es > /ii = Es — Az/2 > Es — Az > ^2 = Es — "iAzjl. Under these 
conditions the dot can be opened and closed for a sequential tunneling current by a spin flip 
induced by an ESR fieldi-W An electron at chemical potential /ii in the left lead and a dot 
electron in state |t) do not have sufficient energy to form the singlet state. If due to an ESR 
induced excitation the dot state is ||), however, less energy is required and an electron in 
lead one can tunnel onto the dot to form the singlet state. Tunneling onto the dot from lead 
two is suppressed by several orders of magnitude if the thermal energy is much lower than 
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the energy gap even if the dot electron is in the excited state ||). At higher temperatures, 
to assure that the singlet can only be formed with the excited dot electron, one can choose 
spin-polarized leads. This may be achieved with several methods, see Ref. |lO| and references 
therein. 

Thus, within these constraints we see that the current in lead one is proportional to the 
probability for the dot being in the excited state, i.e. //(t) oc while the current in lead 
two is proportional to the probability of the dot being in the singlet state, llit) oc ps, see 
Ref.liol 



A. Master equation 

The traditional description of the dynamics of the dot is based on a master equation for 
the reduced density operator of the dot, obtained from the total density matrix ptot(t) by 
tracing over the degrees of freedom of the leads: pdot = T'rieadslPtot]- As usual, we denote 
matrix elements with pnm = I Pdot I (or p„ = (7^|pdot|^)) and include only the three 
relevant dot states n,m E {^,l,S}. We assume the dot and the leads to be uncorrelated 
initially, ptot(O) = Pdot(O) ® Picads(O). Starting from the von Neumann equation for the full 
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density operator ptot = — i [-f^tot; Ptot], the master equation for pdot was derived in Ref. 
using standard methods within the Markov approximation. Further, we will allow for an 
arbitrary (fixed) phase (f of the ESR field which will play an important role in determining 
the spin state. 

In order to eliminate the explicit time dependence emerging from the ESR field, we here 
base our analysis on the dot state in a rotating frame, 

Pdot(t)=e^^*l^><ilpdot(t)e-^-*l^><^l. (2) 

In fact, with the exception of p^ = e*'^*PiT' Pis = 6*'^*Pi5 and the corresponding transposed 
expressions, this transformation leaves almost all matrix elements untouched. 

Along the lines of the derivation in Ref. one finds for the dot state in the rotating 
frame Q a master equation of Lindblad form.«^ It can be written as 

dtPdot = ^Pdot = —^[H, Pdot] + - y^X\LnmPdot, Ll^m\ + [^nm, PdotL\m\) (3) 

nm 

with the time-independent Hamiltonian (in rotating wave approximation) 

if = (A. - c.)|i)(i I + Es\S){S\ - ^(e-nT)(i I + e'1i)(T I) (4) 



and the operators Lnm = y/Wmn\n){m\ describing incoherent transitions between levels m 
and n with a rate Wnm- 

In particular, the four operators Ls-\, Lsi, L^s and L^s describe transitions from and to 
the singlet state and hence correspond to the tunelling of an electron off or onto the dot. 
These four contributions give rise to the current in the leads and are derived from the 
underlying Hamiltonian dTj). The rates are Wsi = Wg^ + Wg^ with W^^ = •yj fi{Es - A^) 
and Wis = Wis + Wfs with Wis = 7/(1 - fi{Es - A,)) where fi{E) = [l + e(-^-^')/*--^] is 
the Fermi function of lead /. Analogously we define the rates Ws-\, W^s, Wsp W^s with 7/ 
and fi{Es). Here 7/ = 27rz//|T/|2 and 7/ = 27rz//| Ti I are the transition rates with density 
of states uj'^ and tunneling amplitude t/"''.«12. In the limit kT <^ A^ we have Wsi ~ 7i and 
Wis ~ 72 5 which resembles the sequential tunneling from lead one onto the dot and into lead 
two. Furthermore, we have Wsi ~ ~ 7i + 72? because we choose fii,fi2 < Es- 

Throughout this paper we assume equal rates for both leads 7i = 72 = 7^ and 71 = 72 = 7'''- 
Finally, we set 7 = 7T = jf ^]^q leads are not spin-polarized, and 7 = 7^^, 7-'' = in the 
case of spin-polarization. 

We simulate stochastically all processes that could be observed in principle, but even- 
tually extract the desired information from those quantities that correspond to the specific 
measurement scheme chosen. Quantum transitions between the dot states may be observed 
by monitoring the current through the dot, which is the starting point for our quantum 
trajectory analysis of the following Sections. 

By contrast, mechanisms for incoherent spin flips (described by the operators and 
Lii) and dephasing mechanisms (described by the projectors L„ = L„„ = \/Wn\n){n\) are 
introduced on phenomenological grounds and not contained in the Hamiltonian (0). The 
(phenomenological) spin flip rates are assumed to satisfy the condition of detailed balance: 
^Ti/^-lT ~ q^^/'^bT^ fjj^g rates Wn are phenomenological dephasing rates: the effect of an 
operator L„„ in (jS)) is to kill coherences between state \n) and the remaining states (at a 
rate Wn), while leaving probabilities unaffected. 

If the coupling to the leads is switched off (by an appropriate choice of the chemical poten- 
tials), the dynamics as described by the master equation ^ is that of a standard decaying 
two-state spin-system. Then the corresponding (intrinsic) relaxation and decoherence rates 
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turn out to be 



1/Ti = W^^ + Wu (5) 
I/T2 = (1/Ti + W^ + W{)/2. 

let us now turn to the dot dynamics: in terms of its coefficients, the time evolution of 
the dot state pdot given by the master equation Q reads 





-^Im(e ''^pit) - (W^ii + Ws^)p^- + W^pi + W^sps, 


(6) 


Pi = 


^Im(e-''^Pll) - {W^ + Wsi)pi + W^uPt + W^sPs, 


(7) 


Ps = 


-{W^s + Wis)ps + W^stPt + WsiPi, 


(8) 


= 


-(i(A, + Vlt)PiT + i^e''^(p^ - pj. 


(9) 


Ps^ = 


-{iEs + Vs^)ps^-i^e'-^Psu 


(10) 


Psi = 


-{i{Es -A, + uj) + Vsi)psi - i^e-^'^psT, 


(11) 



with the effective rates 

Vi^ = ^(VTiT + Wu + Ws^ + Wsi + Wi + W^) = ^{Ws^ + Wsi) + ^, (12) 

Vs^ = liWi^ + Ws^ + W^s + W^s + Ws + W^), (13) 

Vsi = liWn + Wsi + Wis + Wis + Ws + Wi). (14) 

Note that Eqs. (fTIH ITT|) are decoupled from Eqs. (jElini), and the latter are the only ones of 
relevance to us. They enable us to determine easily counting statistics of tunneling electrons 
numerically by means of the quantum trajectory method which we describe in the following 
Sections. 



III. QUANTUM TRAJECTORIES 

One major motivation behind the development of quantum trajectory methods were ex- 
periments with with single quanta. Before these developments, naturally, ensemble experi- 
ments required simple ensemble theories. Matters changed with experiments involving single 
atoms, electrons or ions in traps. Continuously monitoring those systems, single quantum 
jumps became visible to the bare eye. A theory of continuous quantum measurement taking 
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into account continuous measurement records of the observed environment to update the 
quantum state accordingly, were developed, mainly with an eye on applications in quantum 
optics. 

Experiments on the single quantum level have reached solid state devices, as for instance 
electrons in quantum dots. Accordingly, the dynamics of such nano-scale quantum systems 
may be described adequately by quantum trajectories. In fact, it may well turn out that 
these methods are even more useful in solid state devices since the sensitivity of electron 
detectors is typically far better than that for photon detectors, on the single quantum level. 

As will be explained in the following, a quantum trajectory pc{t) describes a subensemble 
of the full (ensemble) density operator p{t), conditioned on a certain (stochastic) measure- 
ment record, here detection events at certain times. In this approach wc determine the 
dynamics of an electron on a quantum dot, conditioned on the measured (stochastic) tun- 
neling current through the dot. 

Quantum trajectory methods have changed remarkably the way we think about open 
quantum system dynamics. While traditionally an open quantum system is described by 
its density operator p{t) as in the last section, quantum trajectories describe open system 
dynamics taking into account certain continuous, stochastic measurement outcomes. In 
other words, with quantum trajectories one determines a conditioned density operator Pc{t), 
reflecting knowledge obtained from a continuous monitoring of the environment. Sampling 
over all these possible measurement records, in other words, ignoring the state of the envi- 
ronment, one recovers the usual full ensemble p{t). We write p{t) — M.[pc{t)] where . .] 
denotes the ensemble mean over all possible measurement records with corresponding prob- 
ability (see below). 

The principle idea is to monitor the environment rather than ignoring, i.e. tracing over it. 
In quantum optics one tries to detect photons emitted from the quantum system of interest, 
here we detect electrons in the leads coupled to the quantum dot. 

In order to illustrate this approach, we consider a simplified open quantum system - 
the generalization to the quantum dot case will be obvious. This model system consists of 
two levels and is coupled to a continuum of states. Excitation is done by some additional 
mechanism, included in the Hamiltonian of the system H. We start with a master equation 
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of type 0, and in this model with a single Lindblad operator L, 

p = Cp = -i[H,p] + ^([Lp,Lt] + [L,pLt]). (15) 

In the following we abbreviate the right hand side of the equation with the superoperator 
Cp. For concreteness, consider L to describe a spontaneous transition from level |1) to level 
|0) with rate W , i.e. L = \/W\{i){l\. We introduce the superoperator S such that 

Kff = {C- S)p = -i[H, p] - ^{L^Lp + pL^L) 
Sp = LpLl 

The latter is referred to as the jump operator since its describes an emission process ac- 
companied by the replacement of the density operator p with the ground state: LpL^ = 
W{l\p\l) |0)(0|. With S such defined, one obtains the quantum jump representatior^ of the 
solution of (fTH|) in the form 

°° ft ftra rt2 

pit) = Y. dtm-1 ■■■ dti (16) 

m=0-^0 Jo Jo 



Pc{t) 

Clearly, the solution p{t) is a sum (or integral, respectively) over any number m of emission 
processes (number of projections onto |0)(0| due to the application of the jump operator 
S), appearing at any times ti,t2, ■ ■ ■ ,tm between zero and the current time t. One has to 
integrate over all corresponding (unnormalized) density operators Pc{t), as apparent from 
expression (fT^ . Thus, one particular quantum trajectory is the normalized density operator 
Pc(t) = pc(t)/tr{pc(t)} which describes the time evolution of the quantum system conditioned 
on the particular measurement record, i.e. conditioned on the number and times of emission 
processes. The quantum trajectory pc(t) occurs with probability tr{pc(t)}. 

The normalized quantum trajectory Pc{t) may be determined directly through the fol- 
lowing prescription: at time t + At the new density operator Pf.{t + At) is obtained in one 
of two ways: 

First, the probability Pjump, to undergo a quantum jump, i.e. to emit a quan- 
tum during the time interval At is equal to the jump rate times the length 
of the time interval times the probability to be in the excited state: -Pjump = 
W{l\pc(t)\l)At =ii{L^Lpc{t))At =ii {Spc{t)} At. If a quantum is emitted (and thus de- 
tected), i.e. a jump has occurred, the conditioned quantum state is the ground state: 



Pc(t + At) = Pjump = |0)(0| = Spc(t)/tr{Spc(t)}. If however, no jump occurs, the new 
density operator is given by 

p.it + At) = p„ = ^^^^n..^;^^^t)} ^^^^ 

as is apparent from the representation (fTBj) . In practice therefore, a quantum trajectory is 
obtained by determining a random number r between zero and one in each time step At: if 
r < Pjump, we set pc(t + At) = pjump, if, however r > Pjump, we set pdt + At) = Pnojump- The 
full ensemble of possible states is thus given by pit + At) = PjumpPjump + (1 — Pjump) Pnojump 
and indeed, one may easily verify that the right hand side equals Cp At as expected from 
the master equation (fT3j) for the full ensemble. 

This branching may occur at any time step and a thus huge ensemble of different quantum 
trajectories may be obtained. As mentioned before, the usual reduced density operator is 
obtained by taking the ensemble mean. In order to obtain counting statistics as in the 
following sections, we simply average over many runs and obtain numerically a distribution 
of jump times as in a real experiment involving a single quantum system. 



IV. COUNTING STATISTICS AND STATE TOMOGRAPHY 

An electron spin on a quantum dot has been found useful as a memory device or a qubit 
for quantum information processing. Readout of the spin state through a tunneling current 
was investigated using a rather restricted parameter regime for which analytical results were 
obtained in Ref. [l^ 

First we want to show, how the analytical results emerge very easily and directly from the 
quantum trajectory approach. Now we consider a regime where we can neglect spin flips, i.e. 
W-^i = = = 0. As in Ref. Q we choose spin-polarized leads, 7-'^ = = Ws-\ = W-^s, 
and 7^ = W. In the limit kT ^ A^ we then have Wsi = W^s = W . The initial state is ||) 
and since no spin flip occurs on the time scale of interest the only processes that happen are 
transitions between ||) and 15). The quantum jump representation (fTH|l of this particular 
solution then reads 



pit) = E / dtrn dtm-1 " " " / dt^ (18) 
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with i,j =l,S and Sijp = W\i){j\p\j){i\. In the regime chosen we can write 

gWeff(tfc-tfc-i)p^(^^__^) = e"^(*''-"*'=-i)pc(tfc-i) and then get 

°° ft Hra 



PW = e-^E/^^-/ dtm-i--- dt^SijSjr--Ssip{0) (19) 
m=o-^o ^0 Jo 

= e-^*f:^K)01P(0)b)(^|. (20) 



Here every operator Ssi describes an electron tunnehng onto the dot from lead one and Sis 
represents the hopping of an electron from the dot into lead two. Since the initial state is 
II), the first transition is ||) \S) and with a second transition back to ||) the first electron 
accumulates in lead two. A third transition to the singlet does not change the number 
of electrons in lead two. For a particular q (number of electrons in lead two) we have to 
consider m = 2q and m = 2g + 1 and the (unnormalized) density operator for a certain q at 
time t is 

P(«.«) = «-""E^K)(4 (21) 

m=2q 

Therefore, the probability to find exactly q electrons in lead two at time t is 

P(,.«)^T.W,*)}^e-2f!(l.J!iL). (22) 

confirming the findings in Ref. Q. 

With the general quantum jump representation p6|). we can overcome the limitations of 
the analytical result, considering arbitrary regimes and investigating the dynamics numeri- 
cally. So far, the proposed measurement scheme allows one to deduce the probability to be 
in either of the two spin states from the current through the dot. A relative phase between 
If) and II), however, cannot be detected. In order to measure the full spin state, there- 
fore, a tomographical measurement setup is required. Here, the freedom to apply the ESR 
field comes into play. We show that while applying an ESR field, phase-sensitive counting 
statistics result, leading to clear identification of the qubit state on the Bloch sphere. As in 
quantum optical setups, the full state could also be obtained with appropriate 7r/2-pulses, 
that effectively change the measurement axis. In this way, not only the (cr^)-component as 
in the original proposal, but also {(Tx) and {(Ty) and thus the full p can be measured. A 
simpler concept, not involving these precise pulses, is to measure the spin state via counting 
statistics of a current through the dot in conjunction with a constant ESR field as we will 
show in the following. For this scheme to be successful it is crucial to control the interaction 
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between dot and leads. We are not interested in the asymptotic, stationary distribution, 
but in the typical time between switching the coupling on and the first (or second, or third 
and so on) electron appearing in lead two. Also, it is not necessary to be able to measure 
the electrons in lead two with a high temporal resolution: one can switch off the coupling 
between dot and lead two after a certain time t and has any time thereafter to collect the 
electrons in lead two. We note that different measurement schemes are possible. Since an 
electron tunneling onto the dot already carries the information about the spin of the dot 
electron, one could abandon lead two altogether and try to monitor the number of electrons 
on the quantum dot, e.g. with a quantum point contact. Our proposal for quantum state 
tomography could be transferred to other setups as well, as for recent experiments.— 

We assume that the dot is in a given initial state at t = 0, when the coupling to the leads 
is switched on. Then we measure the number of electrons tunneling into lead two. According 
to the quantum trajectory approach we calculate the evolution of the density matrix. Every 
jump from 15) to ||) or If) indicates that an electron tunneled out of the dot. At very low 
temperatures, as assumed throughout this paper, the probability of tunneling into lead two 
is close to unity, while tunneling into lead one is very unlikely. 

A single run of the stochastic evolution will display emission processes, i.e. contributions 
to the current, at certain random times. Counting the corresponding number of quanta 
in lead two as a function of time for a large ensemble of quantum trajectories allows us 
to determine the probability P{q,t) of finding exactly q electrons in lead two at time t for 
a given initial state of the dot. Such counting distributions are displayed in the following 
Figures. We stress again that experimentally, it is not required to be able to measure these 
arrival times with a high resolution. One simply switches off the coupling between dot and 
lead two after a given time t. Then one has any time to determine the number of electrons 
in lead two. Our numerical procedure can be applied to any parameter values and any 
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time dependence of the driving ESR field. For the regime chosen in Ref. 
analytical results ()22|1 to a high degree of precision, as will be shown below. 

Let us now turn to the probability distributions P{q,t) of finding exactly q electrons in 
lead two at time t for a given initial state p. As we will show, by employing the ESR field, 
the counting statistics allows to clearly identify the full two-level state, including the relative 
phase. As usual, we choose to parametrise the latter through the coordinates on the Bloch 
sphere: the spin up state ||) corresponds to the north pole with r = 1,9 = 0°, while the 
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FIG. 2: Counting statistics P{q = for the first electron with initial spin up state (dashed 
,r = 1,0 = 0), spin down state (dotted, r = 1,0 = 180°) and the totally mixed state (full 
curve, r = 0). The respective coordinates refer to the Bloch sphere. The inset shows the same 
curves with the counting statistics of the fully mixed state subtracted. Parameters chosen are 
= 27 = 4(A^ - = 2 X 106s-^ T = 2{)mK, = 12T 

spin down state ||) has coordinates r = 1,9 = 180°. The full mixture Po = | (|T)(T| + 
corresponds to the center of the Bloch sphere, r = while coherent superpositions = (|t 
) + e*'^|i))/V2 live on the equator with r = 1,9 = 90°, 0. 

In order to be able to use these counting statistics as a method for spin state tomography, 
the right choice of parameters is crucial. From Eqs. ()fi|7|8|9j) it is obvious that coherences in 
the two-level state can only be transferred to measurable probabilities through the coupling 
introduced by the ESR field of magnitude A^,. On the other hand, a large value of A^ leads 
to Rabi oscillations and thus prevents us to distinguish clearly the two fundamental |t) and 
II) states on a time scale large compared with the Rabi frequency A^/2. Closer inspection 
of Eqs. (jni 13 El) and numerical evidence shows that a good phase sensitivity with preserved 
distinguishability of |t) and |i) is achieved through the choices 

^ ^ Wsi (23) 

Physically, the first condition (on the ESR field strength) means that the spin should not be 
flipped to fast (compared with the measurement time scale Wsi) but still, the ESR field had 
time enough to make the coherences felt. The second condition (on the ESR field frequency) 
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FIG. 3: Counting statistics P{q = for eight coherent superpositions "0 = (IT) + e*'^li))/V2 
along the equator of the Bloch sphere (r = 1, = 90°, various angles 0) and the fully mixed state 
(r = 0, full line). The inset shows the same curves with P{q = of the full mixture subtracted. 
Same parameters as in Fig. |21 

ensures that the method is sensitive to all values of the phase angle 0. 

In Fig. 121 we show counting statistics for the first electron P{q = l,t) to appear in lead 
two. We choose the transition rate 7 = 10^s~^, an experimentally accessible magnetic field 
strength^ of the ESR field ^ 5.16(7, a slightly detuned ESR field frequency — uj = 
5 X 10^s~^, a temperature T = 20mK, and a static magnetic field of strength = 12T. 
For the ESR field to start at zero we choose the fixed phase ip = 37r/2. Furthermore, we 
assume Ti = 10~^s and T2 = 10~^s for the intrinsic relaxation and decoherence times. All 
Figures are calculated with an ensemble of 50000 trajectories. 

The spin down state only allows for electrons to tunnel through the dot, which is clearly 
visible in the counting statistics: if the spin starts off in the spin up state (dotted curve) 
the time to measure the first electron is delayed compared to the mixture and even more so 
compared to the spin down state. Eventually, however, due to the presence of the ESR field, 
a sufficient spin down component will be established allowing electrons to tunnel through 
the dot. Still, both states are clearly distinguishable through their counting statistics. 

Not only are counting statistics useful to distinguish between up and down state. The 
arrival time distribution also differentiates between coherent superpositions and mixtures. In 
conjunction with the ESR field one may even determine the phase of coherent superpositions 
of type '?/' = (II) +e~*'^ ID) /-\/2 as displayed in Fig. 121 The full line corresponds to P(g = l,t) 
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of a fully mixed initial state (r = 0), the dashed and dotted lines correspond to eight coherent 
superpositions along the equator of the Bloch sphere. Clearly, P{q = 1, t) shows different 
behavior for different angles and may thus be used to fully identify the initial state. 

As we have seen, with these choices for the ESR field, not only we can distinguish ||) 
from ID through counting statistics as in Fig. |21 We are in a position to fully determine the 
two-level state - in particular, it is possible to clearly distinguish a coherent superposition 
of It) from II) from the mixture of the two, as shown in Fig. |21 

The insets of Figs. |2 and El reveal an interesting structure underlying the shapes of 
P{q = Once the counting distribution of the full mixture (r = 0) is subtracted, 

statistics of states corresponding to opposite points on the Bloch sphere appear as mirror- 
images of each other, as highlighted in Figs. |3] and El In these Figures we display the counting 
statistics P{q = l,t) for four pairs of opposite initial states along the equator of the Bloch 
sphere and clearly confirm the observations just mentioned. A linear combination of initial 
states leads to a linear combination of counting statistics in the ensemble and thus to this 
symmetry. Still, each curve in itself seems complicated enough to underline the importance 
of our numerical approach. Using the quantum trajectory method, any time dependence of 
the fields and any choice of parameters is possible. 

The more mixed the initial state, i.e. the smaller r < 1 on the Bloch sphere, the closer 
the curve to the curve of the fully mixed state. It is also worth noting that we keep the 
initial phase if of the ESR pulse fixed for all calculations. An average over all possible phases 
would indeed lead to the graph of the fully mixed state, irrespective of the phase (p of the 
initial quantum state. 

A. Higher order statistics and q = 

We close this section by pointing out that also higher order counting statistics (g = 
2,3,4,5) display state-sensitive behaviour - if only less pronounced. This is quite obvious 
since a delayed first tunneling event shifts the starting time for the following electrons. 
For q = 0, the difference between the counting statistics for various initial states is well 
pronounced. In this latter case, however, the curves do not cross which diminishes the 
distinguishability of states and the best choice for that is g = 1. As displayed in Fig. 
ini higher-order counting statistics P{q,t) still distinguishes between the fully mixed state 
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FIG. 4: Graphs taken from Fig. |21for four pairs of opposite states along the equator of the Bloch 
sphere. 




02468 10 02468 10 

7^ 7i 



FIG. 5: Same as Fig. |1] with P{q = l,t) of the fuh mixture subtracted. We clearly see the 
symmetry of the curves for opposite states on the Bloch sphere. 
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FIG. 6: Counting statistics of exactly zero, one, five electrons tunneling through the dot. Evidently, 
if less pronounced, all counting statistics P{q,t) {q = 0, 1, . . . , 5, . . .) carry information about the 
initial quantum state. 

(r = 0) and a coherent superposition {r = 1,6 = 90°). 
B. The role of spin-polarized leads 

The original proposal for the spin state readout was based on spin-polarized leads in order 
to clearly distinguish the two states |t), ||) by a single shot measurement. As the counting 
statistics require an ensemble measurement, our results suggest that spin polarization is not 
required for those - not even advantageous, in fact. In Fig. [7|we display counting statistics 
P{q = l,t) for spin-polarized leads (only spin-up electrons in the leads, i.e. 7-^ = 0). We 
notice only marginal differences compared to the case of unpolarized leads (Fig. OI). 

For large times, it is more likely to observe precisely one electron in the case of unpolarized 
leads. The reason for this behaviour is the fact that for unpolarized leads, there is also the 
possibility that the spin-up electron on the dot (rather than the spin-down electron entering 
the dot) may tunnel out of the dot. Then the dot is in the ground state and therefore closed 
for the tunneling of another electron. It is only after the ESR field had time to populate 
the excited state that a second electron may tunnel through the dot. In fact, it turns out 
that this mechanism is the preferred tunneling event: for the parameters of Fig. |21 we find 
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FIG. 7: Same graph as Fig. 121 Here, however using spin-polarized leads with = 1 x lO^s^^, 
7^ = 0. 

Wis = ly< lOS-i and W^s = 2 x 10^s-\ 

V. RELAXATION AND DECOHERENCE TIMES 

Our proposed setup including the ESR field may be used to determine the intrinsic relax- 
ation time Ti and decoherence time T2 of the qubit in the time domain. Tuning the tunneling 
rate over a wide range (and adjusting the ESR field strength and frequency according to 
conditions (j23p . one can easily see the effect of decoherence and relaxation. In the series of 
graphs in Fig. |Hlwe show counting statistics for the ||) and |t)-state, for the full mixture 
and for two coherent superpositions (states on the equator of the Bloch sphere). Clearly, for 
large tunneling rate (left graph, (a)), all states may be distinguished. The third graph (c) 
shows a regime where decoherence has fully set in: while the states ||) and |t)-state remain 
essentially unaffected, the counting statistics of the coherent superpositions collapse onto 
the curve of the full mixture. In other words, while no relaxation has set in yet, coherences 
between the states ||) and ||) have disappeared. Decreasing the tunneling rate even further, 
counting statistics finally reveal the relaxation time: eventually, the initial states ||) and |t) 
may no longer be distinguished, i.e. the relaxation has taken place. 
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FIG. 8: Measurement of a variety of states (||), ||), coherent superposition) facilitates an estimate 
of T2 and Ti. 

VI. SINGLE SHOT READOUT 

Counting statistics was used in Ref. Q to determine the measurement time and measure- 
ment efficiency. It was shown that after about ten times the tunneling time, the spin state 
on the dot could be determined to be either ||) or ||) with close to 100% efficiency, even if 
only a single measurement is made (provided it was either ||) or ||)). Crucially, these results 
are based on spin-polarized leads. Without spin-polarization the determination of the dot 
state with a single measurement appears problematic. The reason is that an electron may 
tunnel from the quantum dot back into lead one rather than into lead two. In fact, since in 
the case of unpolarized leads there are three ways for an electron to leave the dot: it may 
tunnel from lead one onto the dot and further into lead two (transition \S) or the 

electrons interchange the roles and the one residing on the dot tunnels into either of the two 
leads (15) |t)). The latter two possibilities are almost of same probability and therefore 
with a probability of about one third a tunneling process has taken place without having 
been observed in lead two. In other words, with a probability of two thirds only we can 
claim that without detecting an electron in lead two after sufficient time the spin state on 
the dot was |t). 

In order to overcome this problem we propose to measure the number of electrons on the 
dot with a quantum point contact and recent experimentsi^ show that such a concept may 
work. If an electron tunnels onto the dot it confirms that the state was ||) and if it stays 
there sufficiently long the QPC as electrometer recognizes the charge. Since the projection 
takes place if an electron tunnels onto the dot, the second lead could even be omitted and 
then the dot electron would tunnel into lead one. For clarity in the following argument, 
however, we do not change our setup and leave lead two as before. 
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FIG. 9: With spin-polarized (s-p) leads, a detection of an electron in lead two after about 10 
times the tunneling time ensures to almost 100% that the initial dot state was | |). Without 
spin polarization, only a probability of about two thirds may be achieved due to the possibility 
of the electron to tunnel into lead one. The numerical curve accounts for spin relaxation and 
therefore differs slightly from the analytical result. Without spin relaxation no difference can be 
seen. Parameters as in Fig. |21 respectively Fig. [3 but without ESR- field. 

With the help of the QPC it is possible to measure the state of a single electron spin 
with an efficiency of almost 100% after about ten times the tunneling time, as for polarized 
leads as shown in Fig. ^1 Note that for this scheme to work, it is important that the time 
resolution of the point contact measurement has to be better than the inverse tunneling rate 
to ensure the detection of the tunneled electron on the dot. 

VII. CONCLUSION 

We use quantum trajectory methods to investigate counting statistics of electrons tunnel- 
ing through a quantum dot. We show how an additional ESR field may actively be used to 
perform a full "state tomography". Applying the field during the measurement allows one 
to clearly identify the coherences between the two superposed states. We discuss the rele- 
vance of our findings for determining intrinsic relaxation and decoherence times of electron 
spin states in quantum dots - in the time domain. Based on a quantum point contact we 
propose a scheme for single-shot readout without the need for spin-polarized leads. Similar- 
ities of the investigated quantum dot to three-level systems in quantum optics (e.g. so-called 
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FIG. 10: Measuring the probability to find two electrons on the dot one may determine the dot 
state with 100% efficiency. Clearly, the time resolution has to be better than the typical time the 
two electrons stay on the dot, i.e. better than 7~^. Here P^t'*' the probability that no tunneling 
event occurred in the whole interval [0,t] and Pdot,po('?' ^) probability to find exactly q 

electrons on the dot at time t. The initial state was always ||), parameters as in Fig. |UJ 

"V"-systems) are evident. We underline these connections by applying the quantum jump 
method in order to unravel the dynamics of the full density operator p{t) into subensembles 
pc{t) corresponding to certain measurement records in the leads. Thus, we describe the con- 
ditioned time evolution of the spin state, given a certain measurement record, as in actual 
experiments with single quanta. 

We believe that such connections between methods of quantum optics and mesoscopic 
devices will prove more and more useful in the future as nanotechnology achieves further 
breakthroughs in the coherent manipulation of quantum dynamics in solid state devices. 
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